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Some Cerium compounds exhibit hysteresis cycles with sliarp macroscopic jumps in the magne- 
tization at very low temperatures. This effect is attributed to the formation of clusters in which 
the anisotropy competes with the applied magnetic field. Here, we present a simple model where a 
lattice of ferromagnetically coupled spins is separated in clusters of random sizes and with random 
anisotropy. Within this model, we obtain hysteresis cycles presenting jumps that behave in a similar 
way that the experimental ones, and that disappear when increasing the temperature. The results 
are in good agreement with the hysteresis cycles measured at very low temperatures in CeNii_2,Cua; 
and the comparison with these experimental results allows to discriminate the relative importance 
of the mechanisms driving the thermal evolution of the cycles. 

PACS numbers: 71.27.+a, 75.60.Ej, 61.46.+W, 75.10.-b 



I. INTRODUCTION 

The increasing importance of "inhomogeneous" mag- 
netic systems is due to the fact that inhomogeneities 
(phase separation, local impurities, etc.) are on the basis 
of new emergent magnetic phenomena such as colossal 
magnetoresistance (CMR)ii^ or non conventional super- 
conductivityi^ However, these systems also provide, as it 
is clearly the case of magnetic clustering)^ a direct link 
between the magnetic properties of mesoscopic entities 
and the bulk macroscopical observations. 

Substitutional compounds, in which disorder effects 
are intrinsically present, have been used to tune some 
attractive magnetic effects such as Kondo lattice^ Non 
Fermi Liquid Stated or quantum criticalitiesi^ In par- 
ticular, some well extended theories^ consider systems 
consisting of "magnetic clusters" embedded in a "non 
magnetic" or "paramagnetic" matrix (Griffiths phases), 
to describe the Non Fermi Liquid Stato^. One of the 
more representative examples of such situation is the 
compound CeNii-^^Cu^: For this system it has been 
proposedii a spin cluster-glass state which, by means of 
a progressive percolative process obtained by lowering the 
temperature, finally reaches a long-range ferromagnetic 
state at low temperatures. 

This behavior has been confirmed by two kinds of ex- 
periments. Small Angle Neutron Scattering (SANS) mea- 
surements, which demonstrates the existence of "mag- 
netic clusters" of around 20A, below the freezing temper- 
ature, and low temperature hysteresis cycles, obtained 
in the ferromagnetic state, presenting steps on the mag- 
netization. These steps have been attributed to the 
avalanches of domain flips, as a mesoscopic analogue 
of the Barkhausen noisC j^^i^^ in priority to other possi- 
ble mechanisms proposed in the literature.— i^i^^ In this 
case, it was proved that the "percolated clusters" , reach- 
ing a minimum energy situation, displays the "magnetic 
domain" structure existing in a conventional ferromag- 
net. 

The main characteristics of the observed hysteresis 



loops have recently been presented in an experimental 
articleii and can be summarized as follows. The steps 
are extremely sharp at the lowest measured temperatures 
(~ 100 mK) and strongly dependent on temperature; 
i.e., at 300 mK the steps disappear in the CeNio.6Cuo.4 
compound. The steps appear after the sign inversion 
of the magnetic field and the number of jumps is the 
same in both branches of the hysteresis cycle (increasing 
or decreasing the magnetic field). In Fig. [1] we present 
the hysteresis cycles for two characteristic ferromagnetic 
compounds (CeNio.sCuo.s and CeNio.6Cuo.4). The fer- 
romagnetic long-range order state has been confirmed 
by neutron diffraction at these low temperatures j ^^'^^ At 
this point, we also have to remind that the compositional 
evolution along the series tunes different physical interac- 
tions, such as the Kondo one, which gets enhanced when 
approaching CeNi and also, in particular, the RKKY in- 
teractions, which evolve from the antiferromagnetic be- 
havior of CeCu to a ferromagnetic one when decreasing 
the Cu content This experimental result supports the 
competition of positive and negative magnetic interac- 
tions for the intermediate compositions and, thus, this 
competition has to be accounted for in the model subse- 
quently developed. 

The hypothesis used to explain the emergence of jumps 
in the low temperature hysteresis cycles in this system is 
that magnetic domains developed from a percolative pro- 
cess of static ferromagnetic clusters, reaching a minimum 
energy state at low enough temperatures^^ This mecha- 
nism increases the magnetic correlation length up to val- 
ues that can be detected by neutron diffraction 10'^ A). 
The process is driven by the increasing importance of the 
RKKY interaction as the temperature decreases. This in- 
teraction, then, competes with the local anisotropy, and 
gives rise to a structure of magnetic domains that dis- 
plays an "asperomagnetic" mesoscopic state such as the 
one reported by Coey^ii in the case of amorphous systems. 
The present situation is clearly reminiscent of that case, 
but occurring in crystalline samples. The initial condi- 
tions for such behavior are: random anisotropy for the 
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clusters, disorder and competing magnetic interactions 
between and inside the clusters. A large number of mag- 
netic compounds present special features dealing with in- 
homogeneous states, associated to spin or cluster-glasses, 
phase separations, etd^ i^"*^'^^ The incidence of these fea- 
tures on their magnetic behavior is of great importance 
in determining and characterizing new ground states of 
matter (Quantum Critical Points, etc.) at low tempera- 
tures. 

A first attempt to theoretically describe the magneti- 
zation jumps in the hysteresis cycles of the CeNii_3;Cua; 
alloys was previously developed by using a spin Hamilto- 
nian, including a spin-spin interaction term, an applied 
external field and a random local ficldj^ This random 
local field was intended to simulate the effect of a local 
anisotropy, but instead of considering just an axis of easy 
magnetization, there was also a preferred sense, what 
means that the anisotropy energy was minimum when the 
spin was aligned with the local random field, and maxi- 
mum when pointing in the opposite sense. The hysteresis 
cycles so obtained presented a number of jumps strongly 
dependent on the value of the spin, (assumed equal or 
bigger than one half), and on the number of possible dis- 
crete orientations for the random field. Even if this model 
describes the qualitative features of the phenomena, the 
jumps in the hysteresis cycle are equally spaced, as both 
the spin and the anisotropy can only have discrete val- 
ues. Moreover, in that case, the anisotropy was locally 
described at each site. Thus, the model is not well suited 
to represent a system of clusters. 

We propose here an alternative model that seems to be 
better adapted to the description of the studied effects. 
On one hand, this model is simpler than the previous 
one, being done that now the spins are described by an 
Ising model and have just two possible orientations. On 
the other hand, anisotropy can now have random orienta- 
tion, and clusters of random different sizes with different 
anisotropy directions are considered. Then different clus- 
ters have different anisotropy orientation, but inside each 
cluster all spins are subjected to the same anisotropy di- 
rection, which is not related to the cluster size. This way 
of describing the anisotropy field directly arises from the 
model proposed to describe the magnetic behavior of the 
CeNii_2;Cu2: compoundsiii and is a useful tool to de- 
scribe bulk properties arising from a microscopic dynam- 
ics. Finally, the temperature effects are not restricted to 
the thermal activation of the spins, but a phenomenolog- 
ical description of the percolative process is also consid- 
ered. We also remark as a difference with the previous 
model that each cluster is now characterized by a true 
anisotropic term; i.e., by an axis of easy magnetization, 
whose direction is random, and not by a random field, as 
in ref. — . We should recall at this point that, although 
the random anisotropy model developed in the present 
work provides a description closer to the empirical point 
of view, there exists much further literature on random 
field models, including extensive reviewsj ^^'^^ However, it 
is also true that both kind of models are closely related^^ 
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FIG. 1; Hysteresis cycles of two representative compositions 
(CeNio.sCuo.s and CeNio.6Cuo.4) measured at very low tem- 
peratures. Note that the jumps take place upon reversal of 
the magnetic field and the number of jumps is the same in 
both branches of the hysteresis loop. 

and, even, both of them belong to the same universality 
class.— 

Therefore, we have undertaken a numerical study to 
model the hysteresis cycles of CeNii_a;Cua; as a function 
of temperature and we present here the results of our 
Montecarlo calculations. In the next section we will de- 
scribe the model and the simulations, while in section III 
wc present and discuss the results. 

II. DESCRIPTION OF THE MODEL 

The model we consider is a three dimensional Ising sys- 
tem on a simple cubic lattice. The lattice is separated 
in clusters of random mesoscopic size, and each cluster is 
characterized by a random anisotropy direction. The in- 
teraction between the spins is mainly ferromagnetic but 
we allow a given concentration of antiferromagnetic in- 
teractions to represent the disorder and the fact that the 
magnetic order is different in the limit of high Ni or high 
Cu concentrations. Then, the Hamiltonian of the system 
reads: 

i.j n i(^n i 

(1) 
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where 5*^ are Ising spins that can take the values ±1, so 
describing a system with spin 1 /2. This situation is repre- 
sentative for a 00"^+ ion in an orthorhombic environment, 
as its J = 5/2 ground state is spht into three doublets 
and, thus, only one doublet is relevant at very low tem- 
peratures. The spins interact among them through the 
magnetic interactions Jij = ± J, the plus sign represent- 
ing ferromagnetic interactions and the minus sign antifer- 
romagnetic ones. The value of the exchange constant is 
used in the present work as the scaling factor to present 
the results on other magnitudes; i.e., the anisotropy, the 
temperature and the magnetic field will be measured in 
units of J. We also consider that the lattice is divided 
in clusters of random sizes, and that each cluster, n, is 
characterized by an anisotropy field An whose absolute 
value is assumed to be constant. A, while the orientation 
with respect to the z-axis, 0„, is random. We remark 
at this point that the anisotropy strength and direction 
is assigned to each individual spin, and that the calcu- 
lations are made on the basis of these individual spins. 
Therefore, the energy barrier for spin flip is related only 
to the local strength of the anisotropy, A cos </>„ , the ex- 
change energy, J (times the number of neighbors), and 
the applied field, Ji. So, differently to the previous mod- 
els (including rcf. [23i ) here the system is separated in 
clusters of random sizes, and the second term correspond 
to an anisotropy interaction characterized only by the di- 
rection (and not the sense) of the anisotropy axis. More- 
over, we are trying to describe a disordered compound 
in which there is a strong competition between positive 
and negative magnetic interactions. We must recall here 
that CeCu is antiferromagnetioi^ while a ferromagnetic 
ordered phase builds up at low temperatures for the Ni- 
richer alloys, as a consequence of a percolating process 
of magnetic clusteringJ^iii Hence, as we are interested 
in the hysteresis cycles at very low temperatures, we as- 
sume that in the Hamiltonian (1) the interactions Jy 
arc mainly ferromagnetic, Jy = -I- J, and we introduce a 
small, but finite, concentration of AF links, Jy = — J. 

Then, within this model, we have performed a Monte 
Carlo simulation on a 3-dimensional lattice, at zero and 
finite temperatures, considering a cubic lattice of dimen- 
sion N X N X N with N ~ 80. We considered a number 
of clusters of the order of N, because we have observed 
that a good qualitative description of the experimental 
results is obtained for relatively big clusters and a small 
dispersion in size. The strength of the anisotropy has 
been taken as A = 3J, and the random orientation be- 
tween and TT. Then, clusters are generated according to 
the following recipe: a) an arbitrary number of spherical 
clusters R is assumed, being R of the order of N, b) the 
center of each cluster is determined at random within the 
sample, c) The radius of the cluster is chosen at random 
with uniform probability in a range {rmin,i"max}, d) to 
prevent superposition of clusters, sites are labeled as be- 
longing to a cluster once its anisotropy is determined. 
Eventually some sites may be out of any cluster, but the 
number of such sites not belonging to any cluster is much 
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FIG. 2; Two dimensional section of the lattice showing the 
clusters and the relative angle of the anisotropy, between —tt 
and — TT. For this figure we use Vmin = 0.2A^ and Vmax ~ 0.57V 



smaller than the total number of sites. To illustrate the 
cluster geometry, in Fig. [2] we show a section of a lat- 
tice exhibiting a typical 2D cluster distribution. It is 
evident that some big clusters are present and they are 
the ones that determine the strength of the jumps in the 
hysteresis cycle. In the present simulation we have taken 
R = N and relatively big clusters, 0.25 < < 0.45 

and 0.4 < r,nax < 0.65. 

We have also determined the cluster size distribution. 
As the samples are rather small to perform a statistics, 
we have studied 400 samples, and the results are plotted 
on Fig. [3l The distribution is represented in a semi-log 
plot, and the behavior is almost an exponential law. The 
average size of the clusters is 18 sites but with a relatively 
high variance, of the order of 44. That means that there 
are a big number of small clusters filling the space among 
a few large clusters. Those large clusters are enough to 
cover most of the lattice, and we can see from the figure 
that the frequency of large clusters (size bigger than 10^ 
sites) is of the order of the number of samples. 

We should keep in mind that, as the anisotropy term in 
the Hamiltonian (1) is quadratic, the contribution of this 
term to the total energy is independent of the spin being 
S" = ±1. This fact implies that, if one considers a dynam- 
ics leading to the state of minimum energy, there would 
be no effect of this term. Hence, we assume that, in order 
to change the spin direction in one site, the system has 
to overcome an energetic "barrier"; i.e., it needs to have 
enough energy to outdo the strength of the anisotropy 
field. Therefore, at zero temperature, we consider that 
if the spin is in a state of negative energy, it will re- 
main there independently of any consideration about a 
possible (more negative) energy of the state with oppo- 
site spin. Just when the energy becomes positive, the 
spin will reverse its direction. Finally, as it was stated 
above, a percentage of the links Jij are taken to be anti- 
ferromagnetic. In the results presented in this work, this 
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FIG. 3; Distribution of the size of the clusters (measured as 
the number of sites of the cluter) evaluated over 400 samples, 
with the same parameters of Fig. [5] 

percentage is considered to be 5%. 

Two different processes that influence the temperature 
evolution of the jumps observed in the hysteresis cycles 
will be considered: 

a) Thermal activation: This is the usual process where 
the spin always changes its direction if its energy is posi- 
tive but, even if the energy of the spin Si on the i-site is 
negative, — i?, we assume that there is a finite probability 
ig-B/T the spin changes its orientation. 

b) Cluster percolation: In this case, we will have to 
account for this effect in a qualitative way, as we do not 
have a precise mathematical model on the temperature 
evolution of the percolation paths. However, what we 
know from the experimental resultSfii is that the per- 
colation process yields an increase of the ferromagnetic 
correlation length when lowering the temperature. Thus, 
in our calculation, we can associate the different stages of 
this percolating process with different cluster size distri- 
butions; each of the distributions centered on larger sizes 
as the temperature goes down. 

III. RESULTS AND CONCLUSIONS 

In Fig. m we present a typical hysteresis cycle obtained 
at T = using the model described above. The hystere- 
sis cycles obtained with the simulation exhibit some re- 
markable coincidences with the experimental results (see 
Fig. [T]). Not only the number of jumps is the same when 
inverting the magnetic field but the cycle presents the 
same inversion symmetry with respect to the origin, char- 
acteristic of the experimental results. And those features 
are independent of the number of jumps obtained by us- 
ing different seeds for the random number generator in 
the simulation. Also, the disorder introduced with the 
elimination of some ferromagnetic links produce a curva- 
ture of the magnetization curve for low values of the mag- 



FIG. 4: One example of the hysteresis cycles obtained by 
numerical simulation, with the parameters described in the 
text. The magnetic field is measured in units of J and the 
magnetization is relative to the saturation value. 




FIG. 5: Effect of temperature activation only: Variation of a 
simulated hysteresis cycle as a function of temperature. The 
steps practically disappear for a temperature of the order of 
T ~ 0.5J. The parameters (size of the lattice, anisotropy, 
disorder, size distribution of the clusters) are the same as 
those used for obtaining Fig. 2] but different seeds have been 
taken for the random number generator. 



netic field, W, that qualitative coincides with the exper- 
imental results and was not observed in previous calcu- 
lationsi^ This qualitative agreement between simulation 
and experiments strongly suggest that the existence of 
clusters with different anisotropy directions is at the ori- 
gin of the jumps in the hysteresis cycles, as well as other 
properties of the CeNii_a;Cua; compounds. Each of the 
jumps of the hysteresis cycle corresponds to an avalanche 
process where the spins (of one or more clusters) align 
with the applied magnetic field. The bigger the size of 
the cluster, the bigger the avalanche and the lower the 
number of jumps. So, one can conclude from the simula- 
tions that the clusters are of significant, mesoscopic size, 
as if one reduces too much the size of them, the number 
of jumps should be much larger than in the experimen- 
tal results. We remark that the number of jumps is so 



5 



0.4 



0.2 



m 0.0 



-0.2 



-0.4 



' 1 ' 
CeNio.6Cuo.4 




1 






— T = 100 mK 




~ T = 200 mK 




~ T = 300 mK 




^ T = 500 mK 


1 


1 



-2000 -1000 1000 



H(Oe) 

FIG. 6: Temperature evolution of the left branch of the hys- 
teresis cycle for a representative composition (CeNio.6Cuo.4). 
Note that the coercive field almost does not change whilst the 
jumps are still noticeable, and it clearly decreases for higher 
temperatures, once the jumps have disappeared. 




FIG. 7: Effect of the cluster percolation on the hysteresis 
cycles. The different stages of the percolation process are 
correlated with cluster distributions around different sizes: 
the larger the size, the lower the temperature. The rest of the 
parameters (size of the lattice, anisotropy, disorder) are the 
same as those used for obtaining Fig. [S] 



low that it is not possible to make an statistical analysis, 
and if one makes statistics over different samples the size 
of the jumps is so different among different samples that 
the effect would disappear. However, we have performed 
simulations with different size of the samples and/or the 
clusters, and with different seeds for the random number 
generator and the results are consistent with the above 
mentioned conclusions 

We have also performed finite temperature calcula- 
tions. First we have just consider thermal activation: 
an example of the temperature variation of the hystere- 
sis cycles is presented in Fig. [5l where one can see that 
the jumps in the hysteresis cycle disappear as the tem- 
perature increases, in very good agreement with the ex- 
perimental data, presented in Fig. [51 In the simulation, 
the temperature is measured in units of J. So, compar- 
ing with the experimental results, where the jumps are 
fully wiped out at temperatures of the order of SOOmif, 
we can estimate the value of the exchange interaction as 
being of the order of \K. This is a too small value for the 
exchange interaction, but this result can be attributed to 
the fact that in the simulation the number of spins in each 
cluster is several orders of magnitude smaller than in the 
physical system. In addition, we should keep in mind 
that this system exhibits some characteristics of Kondo 
behavior, so that the value of the magnetic moments at 
very low temperatures are considerably reduced, as ob- 
served by neutron diffractionJ^ii^ 

Then, we have tried to show the effect of the perco- 
lation process. As stated before, the different stages of 
this percolation process can be qualitatively correlated 
with cluster distributions centered around different sizes: 
the lower the temperature, the larger the size of the clus- 
ters. These qualitative results are presented in Fig.[7]and 



the observed behavior shows clearly that the jumps are 
smoothed as the mean radii of the clusters is reduced; 
i.e., the temperature is increased. This result is an indi- 
cation that the cluster percolation process is also able to 
account for the temperature evolution of the experimen- 
tal hysteresis cycles. 

Although ab-initio considerations arising from the pro- 
posed model do not allow to clearly distinguish which one 
of the mechanisms is more relevant in order to define the 
variation with the temperature, the analysis of the re- 
sults concerning both of these mechanisms may give us 
relevant information in this sense. If we look at the coer- 
cive field variation arising from both situations, we will 
realize that it is much larger when considering the ther- 
mal activation process than when taking into account 
the percolation one. Comparing now with the experi- 
mental results presented in Fig. [S] we can observe that 
the coercive field do not hardly change whilst the jumps 
in the hysteresis cycle are still noticeable (T < 300 mK), 
whereas it strongly decreases with the increasing temper- 
ature once the jumps have completely disappeared. This 
fact is a clear indication that the percolative process is 
the dominant one at the lowest temperatures, but once 
a "conventional" hysteresis cycle is reached, the thermal 
activation becomes predominant. 

So, keeping in mind the limitations of the present sim- 
ulation, performed with a relative small number of spins, 
the model proposed here can well account for the strik- 
ing properties of the hysteresis cycles of the CeNii_a;Cua; 
compounds, describing both the jumps of the magnetiza- 
tion at zero temperature and the extinction of the jumps 
when the temperature increases. A number of facts indi- 
cates that the model presented rcfiects in an appropriate 
way the experimental evidence and the proposed scenario 
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of percolative ferromagnetic clusters for this series: 

a) The estimated low value for the exchange energy is 
a signature of the Kondo interaction. 

b) Some amount of disorder in the magnetic interac- 
tions is necessary to account for the curvature of 
the magnetization curves in the low field range. 

c) The clusters must be of mesoscopic size in order 
to obtain a small number of jumps, as observed 
exp er iment ally. 

d) The formation of these clusters of mesoscopic size 
is intrinsically related to the disorder in the interac- 
tions and the small value of the magnetic moments. 

e) The temperature evolution of the hysteresis cycles 
is dominated by the percolation process at the low- 
est temperatures and by the thermal activation at 
higher ones. 



These points are critical for an understanding of the fer- 
romagnetic cluster percolative model. In similar series 
based in other Rare-Earths, such as Nd or Tb, no evi- 
dence of cluster creation has been foundi^ Therefore, we 
can conclude that the calculations presented here provide 
a remarkable support for the model proposed in refer- 
ences [l3 and [iTI. which, in fact, could be a valid descrip- 
tion for a large number of strongly correlated electron 
systems with disorder effects. 
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